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Octacalcium phosphate (OCP) is prov- 
ing to be an important intermediary in 
the formation of tooth and bone mineral 
and various pathological calcifications. 
Before this mineral can form, its solubil- 
ity product must be exceeded. Thus, a 
knowledge of its precise values under 
various conditions is required for a basic 
understanding of calcification processes. 
The methodology suitable for measur- 
ing the solubility of metastable phases 
was developed and used to determine 
the negative logarithms of the solubility 
products of OCP, p^,p(OCP), at 4, 4.8, 
6, 18, 23.5, and 37 °C. This methodology 
includes (1) the use of high solid-to-liq- 
uid ratio, 10 mg/mL, to minimize the 
effects of hydrolysis, (2) frequent sam- 
pling during equilibration to detect 
possible effects of hydrolysis, (3) equili- 
bration from supersaturation and from 
undersaturation, and (4) equilibration in 
the absence and presence of a C02-con- 
taining atmosphere. 



The resulting pX'sp(OCP) values are 
48.3±0.2, 48.3 ±0.2, 48.2, 48.3, 48.4±0.1 
and 48.7±0.2 at 4, 4.8, 6, 18, 23.5, and 
37 °C. A 5.5% CO2 atmosphere did not 
change the apparent pA',p(OCP) value 
significantly. The value of p^sp(OCP) 
obtained by approaching equilibrium 
from supersaturation was essentially the 
same as that from undersaturation. 

The effects of (1) the use of different 
ionic models, (2) OCP hydrolysis, and 
(3) differences in equilibrium constants 
on the apparent p/TspCOCP) values are 
described; the latter two contribute 
significantly to the differences in 
p^sp(OCP). 
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Introduction 



Octacalcium phosphate (OCP), Ca8H2(P04)6- 
5H2O, is frequently encountered in calcium phos- 
phate systems more basic than dicalcium phosphate 
dihydrate (DCPD), CaHP04-2H20. It forms and 
hydrolyzes rapidly under physiological conditions 
[1]. There are many indications that it is the precur- 
sor in the formation of tooth and bone minerals and 
pathological calcium phosphates [2,3]. The struc- 
ture of OCP has been determined [4] and its solu- 
bility reported [5-8]. The recently pubhshed values 
[5,6] for pJ^sp(OCP) (pir,p(OCP)=-log K,^{OCP), 



where ^sp(OCP) is the solubility product of OCP) 
are in variance with older values [7,8] as shown in 
table 1. Thus, it is important to know the precise 
thermodynamic solubility of OCP in the form of 
the solubility product [p.Srsp(OCP)] in order to as- 
sess the degree of saturation of biological fluids 
with respect to OCP. We report here values for 
p^,p(OCP) at 4, 4.8, 6, 18, 23.5 and 37 °C. One of 
the main difficulties in obtaining the solubility of 
OCP is that OCP is not the most stable phase and 
will hydrolyze to other thermodynamically more 
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Table 1. Solubility products of OCP, as p^,p(OCP), at various temperatures, this work and 
literature values. Values in parentheses were recalculated using our ionic equilibrium model 



Temp. 
CO 



This work 



4 


48.3±0.2 


4.8 


48.3±0.2 


6.0 


48.2 


18 


48.3 



Nancollas 
et al. [5,6] 



Madsen [7] 



Moreno 
et al. [8] 



23.5 



25 



37 



45 



48.4±0.1 



48.7±0.2 



49.6±0.2' 
(49.0±0.2)'= 

49.3+0.2" 
(49.1 ±0.3)" 

49.8±0.3* 
(48.6±0.4)'' 



46.97" 



48.46' 
(48.0)" 



" Calculated by Shyu et al. [5] and Heughebaert and Nancollas [6]. 

" Recalculated by using our ionic models and constants. 

' Recalculated by using our ionic constants; when ionic activity coefficients were calculated 

from the extension of the Debye-Hiickel limiting law, the p^=48.97±0.25; from the Davies 

equation pA'=49.03±0.23. 



Stable calcium phosphates. Therefore, procedures 
suitable for measuring the solubility of metastable 
phases were developed and used. The methods 
minimize the effect of OCP hydrolysis and include 
(1) the use of high solid-to-liquid ratio, (2) frequent 
sampling during equilibration to detect possible ef- 
fects of hydrolysis, (3) equilibration from supersat- 
uration and from undersaturation, and (4) 
equilibration in the absence and presence of a CO2- 
containing atmosphere, since carbonate slows 
down the hydrolysis. We have also studied the ef- 
fect of OCP hydrolysis on the apparent solubility 
product and examined the possible causes for the 
differences in OCP solubility constants in the litera- 
ture as compared to the values reported here. 



Materials and Methods 
Materials 

OCP was prepared by two methods: (1) Hydrol- 
ysis of DCPD slurry in distilled water with addi- 
tion of water dropwise over a period of about 3 
months at room temperature [9]; this sample was 
designated as OCP-A and used for most of the 

Study unless specified otherwise. The DCPD that 

was used for preparing the OCP was prepared by 
ammoniating an aqueous solution initially saturated 



with monocalcium phosphate monohydrate 
(MCPM), Ca(H2P04)2-H20, and DCPD, (i.e., at the 
singular point of MCPM and DCPD) according to 
the procedure described by Moreno et al. [10]. The 
OCP had a Ca/P ratio of 1.33±.02 and yielded a 
characteristic OCP x-ray diffraction pattern. (2) 
Growth of OCP by the constant composition 
method at pH=6[ll]. This preparation was desig- 
nated as OCP-B. 

Analytical Methods 

Two spectrophotometric methods were used to 
determine the concentration of calcium and phos- 
phate: calcium as the Ca-Arzenazo III complex 
[12] and by atomic absorption [13], phosphate as 
vanado-molybdate [14] complex or the phospho- 
molybdate malachite green complex [12]. The esti- 
mated errors were ±1.5% of the amounts 
analyzed. The pH was measured with a combina- 
tion glass-calomel reference electrode with esti- 
mated errors of ±0.008 pH units. 

Equilibration Procedures 

The saturated solutions of OCP were obtained 

by allowing OCP crystal dissolution or growth to 
proceed to equilibrium in solutions initially under- 
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saturated or supersaturated with respect to OCP. 
Crystals of OCP were added to various concentra- 
tions of phosphoric acid (10 mg/mL), when equi- 
librium was approached from undersaturation, or 
to solutions containing potassium dihydrogen phos- 
phate and calcium nitrate, when equilibrium was 
approached from supersaturation. The suspensions 
were stirred with a magnetic stirrer, and the exper- 
iments were carried out either in covered jacketed 
Pyrex cells through which thermostated water was 
circulated or in sealed plastic tubes which were 
mounted in a constant temperature bath at 4, 4.8, 6, 
18, 23.5, or 37 °C. In some of the experiments the 
nitrogen or 5.5% C02/94.5% N2 mixture presatu- 
rated with water was bubbled through the solu- 
tions. The equilibrations were allowed to proceed 
for different periods of time (15 min to 3.5 months). 
Samples were withdrawn periodically from the 
suspensions and either centrifuged at 15,000 rpm in 
an Eppendorf centrifuge (No. 5414) or filtered 
through a 0.22 jum filter (Millipore, Bedford, MA).' 

The pH was measured just before taking the 
samples, after separation of the solids, and, in some 
cases, continuously at the initial stage. The calciimi 
and phosphate concentrations of the supematants 
or filtrates were determined as described above. 

The ion activity products for OCP [/P(OCP)] 
were obtained from the calcium and phosphate 
concentrations, pH and ionic strengths using a 
computer program [15] which calculates the ion 
activity coefficients through the use of the ex- 
tended Debye-Hiickel equation [19] or Davies 
equation [20] and takes into account the ion pairs 
CaOH+, CaHP04, CaH2P04+, CaCOj and 
CaHCOj"^ [16], and dissociation of H3PO4 and H2O, 
(Appendix) as follows: 

IP(pCF) = (Ca='+)\H+)(P04'-)^ 

Here parentheses indicate ionic activities. The 
A:sp(OCP)=/P(OCP) when the solution is in equi- 
librium with OCP. Unless specified otherwise, the 
extended Debye-Hiickel equation was used. 



' Certain commercial materials and equipment are identified in 
this paper to specify the experimental procedure. In no instance 
does such identification imply recommendation or endorsement 
by the National Bureau of Standards or the ADA Health Foun- 
dation or that the material or equipment identified is necessarily 
the best available for the purpose. 



Results 

The compositions of solutions equilibrated with 
two different preparations of OCP and p/P(OCP) 
values as a function of time at 37 °C in a nitrogen 
atmosphere are presented in tables 2 and 3 (table 2 
for solutions initially undersaturated, and table 3 
for solutions initially supersaturated with respect to 
OCP). The results are similar for both preparations 
whether or not the initial solutions are under or 
super saturated with respect to OCP, except after 
10 days. Equilibrium was obtained rapidly from 
both undersaturation and supersaturation (within 
15 to 20 min) as indicated by constancy of the val- 
ues of p/P(OCP) and pH of the solutions. In the 
equilibration of OCP-A from imdersaturation, the 
p/P(OCP) was constant for about 10 days (depend- 
ing on the initial phosphoric acid concentration) 
and then increased slowly after that. This increase 
is attributable to the hydrolysis of OCP to a more 
basic product which has lower solubility than that 
of OCP under the conditions studied. The pH of 
the solution increased rapidly in the first 15 minutes 
due to the dissolution of OCP, remained constant 
for 4.5 h, and then decreased slowly over the range 
of time studied, as indicated in table 2. The latter 
decrease can only be due to the hydrolysis of OCP 
to a non-stoichiometric "apatite." The hydrolysis 
produces phosphoric acid approximately as 
follows: 

(10-;<:)Ca8(P04)6H2.5H20^8Ca,o-.(P04)6(OH)2_2, 

-K12-6x)H3P04+(34-f llx)H20. 

The Ca/(P— Po) ratios in solutions should indicate 
the dissolution ratios. These ratios, although vary- 
ing from the ratio for congruent dissolution 
(4/3), decreased with time, indicative of a hydroly- 
sis process which increased the phosphate concen- 
tration in the solution. After equilibration of 
sample OCP-A for 3.5 months, only apatite lines 
could be detected in the x-ray powder diffraction 
patterns, and the p/P(OCP) was 50.1. When equi- 
librium was approached from supersaturation, the 
pH of the solution decreased rapidly from 5.7 to 
5.47 in the first 20 minutes, as indicated in table 3, 
due to the precipitation of calcium and phosphate 
from solution and crystal growth of OCP. After 
that, the pH of the solution decreased slowly, and 
yet the pZP(OCP) remained constant. The slow de- 
crease in the pH of the solution is again indicative 
of OCP hydrolysis to an apatite-like phase. Al- 
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Table 2. Various times of equilibration of two batches of OCP with 0.457 mmol/L HjP04 solutions from undersaturation 
at 37 "C and under N2 atmosphere 







Composition of solution 






Time of 


pH 


[Ca] 




[P] 


Ca 


p/P 


equilib. 




mmol/L 




mmol/L 


P-Po 




OCP-A 





3.35 


0.0 




0.457 






15 min 


6.35 


0.614 




0.994 


1.14 


48.86 


31 min 


6.32 


0.616 




0.993 


1.15 


48.98 


62 min 


6.30 


0.630 




1.032 


1.10 


48.98 


2h 


6.36 


0.644 




1.052 


1.09 


48.68 


4.5 h 


6.35 


0.670 




1.092 


1.06 


48.68 


23 h 


6.16 


0.806 




1.367 


0.89 


48.88 


10 d 


5.43 


2.91 




5.37 


0.59 


48.78 


2 months 


4.75 


7.29 




13.91 




49.69- 
av 48.83 


OCP-B 





3.33 


0.0 




0.457 






15 min 


6.38 


0.685 




0.959 


1.36 


48.61 


33 min 


6.37 


0.703 




0.982 


1.34 


48.58 


65 min 


6.40 


0.694 




0.972 


1.35 


48.65 


2h 


6.41 


0.701 




0.996 


1.30 


48.38 


4.5 h 


6.43 


0.685 




0.969 


1.34 


48.36 


23 h 


6.33 


0.767 




1.070 


1.25 


48.50 


10 d 


6.10 


1.022 




2.912 


0.42 


48.32 


2 months 


4.96 


6.57 




13.52 


0.50 


48.84" 
av 48.49 



" Were not included in solubility calculation. 



Table 3. p/P of OCP from supersaturation at 37 °C 







Composition 


of solution 






Time of 


pH 


[Ca] 


[P] 


Ca— Cao 


pIP 


equilib. 




mmol/L 


mmol/L 


P-Po 







5.70 


4.43 


3.33 




47.59 


20 min 


5.47 


4.41 


3.38 


-0.40 


48.68 


Ih 


5.46 


4.37 


3.36 


-2.00 


48.75 


2h 


5.45 


4.35 


3.39 


-1.33 


48.80 


4h 


5.43 


4.42 


3.43 


-0.10 


48.86 


22 h 


5.41 


4.60 


3.73 


0.43 


48.79 
av 48.78 





5.70 


4.80 


3.28 




47.50 


42 min 


5.35 


5.53 


4.00 


1.01 


48.71 


2d 


5.30 


5.61 


3.98 


1.16 


48.94 


3d 


5.27 


6.02 


4.44 


1.05 


48.87 


5d 


5.22 


6.67 


4.95 


1.12 


48.83 


12 d 


5.07 


8.42 


6.77 


1.04 


48.86 
av 48.84 
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though the hydrolysis of OCP-A seemed to occur 
early, as indicated by the decrease in pH and 
ACa/AP in 23 h, the effects of hydrolysis on the 
value of pIP(OCF) was observed only after 10 
days in the equilibration from undersaturation 
(table 2). The hydrolysis of OCP-B was slow com- 
pared to that of OCP-A, as indicated by the change 
in pH of the solution (table 2); therefore, the varia- 



tion of p/i'(OCP) with time was not as great as in 
the case of OCP-A. 

The resuhs of pK,p(pCF) under a 5.5% COj at- 
mosphere or under air (tables 4 and 5) fall in the 
same ranges as those under N2 (tables 2 and 3). The 
hydrolysis rate, as indicated by decreases in pH and 
ACa/AP in tables 2, 4, and 5, depends on the CO2 
content in the atmosphere and decreases in the or- 



Table 4. p^,p of OCP under 5.5 volume % of CO2/94.5 % Nj atmosphere from undersaturation and at 37 "C 







Composition 


of solution 






Time of 


pH 


[Ca] 


[P] 


Ca 


P^,P 


equilib. 




mmol/L 


mmol/L 


P-Ps 







3.196 





0.509 






1,5 h 


5.996 


1.27 


1.53 


1,24 


48.80 


2.0 h 


5.937 


1,31 


1,57 


1,23 


49.00 


28.0 h 


5.928 


1.30 


L69 


1.10 


48.97 


2d 


5.949 


1.30 


1,73 


1.06 


48.84 


5d 


5.914 


1.50 


1.99 


1.01 


48,61 


7d 


5.930 


1.49 


2,04 


0,97 


48,52 


9d 


5.812 


1.63 


2,J5 


0,99 


48,87 


12 d 


5.764 


1.80 


2,65 


0,84 


48.70 


16 d 


5.675 


1.95 


2.99 


0.79 


48,86 
av 48.80 





2.962 





0.935 






1.5 h 


5.808 


1.42 


2.07 


1.25 


49,14' 


2.0 h 


5.808 


1.52 


2.16 


1.24 


48.99 


28,0 h 


5.814 


1.46 


2.16 


1.19 


49.03 


2d 


5.841 


1.49 


2.22 


1.16 


48.83 


5d 


5.817 


1.64 


2.36 


1,15 


48,72 


7d 


5.853 


1.59 


2.32 


1,15 


48.63 


9d 


5.773 


1.66 


2.35 


1,17 


48,92 


12 d 


5.742 


1.86 


2.84 


0.9S 


48,68 


16 d 


5.695 


1.99 


3.03 


0.95 


48,72 
av 48,82 



' Was not included in solubility calculation. 



Table S. pK^ of OCP under air atmosphere at 37 °C 







Composition 


of solution 






Time of 


pH 


[Ca] 


[P] 


Ca 


pK>, 


equilib. 




mmol/L 


mmol/L 


P-Po 












0.463 






0.5 h 


6.35 


0.606 


1.073 


0.99 


48.80 


Ih 


6.33 


0.621 


1.083 


1.00 


48.84 


2h 


6.32 


0,650 


1.083 


1.05 


48.81 


28 h 


6.10 


1,056 


1.756 


0.82 


48,44 
av 48.72 









0.463 






2h 


6.33 


0,630 


1,139 


0.93 


48.74 


2d 


5.83 


1.384 


2.500 


0.68 


48.86 


3d 


5.68 


1,74 


3.424 


0.59 


48.86 
av 48.82 
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der N2>air>5.5% CO2 due to the inhibitory effect 
of carbonate on hydrolysis. On the other hand 
there is no observable effect of carbonate on the 
p^sp(OCP) presumably because carbonate cannot 
be incorporated into the OCP crystals. 

The solubility products under different initial so- 
lution compositions at 4, 6, 4.8, 23.5, and 37 °C are 
shown in tables 6-11 together with the equilibrium 
times, composition of solutions and ACa/AP, The 
results indicate that hydrolysis is slower when the 
initial solution has lower pH (i.e., higher phospho- 
ric acid concentration), lower temperature, or the 
initial solution is supersaturated. There are indica- 
tions of lack of attainment of saturation in the early 



stages of the equilibration in the presence of CO2 
(table 4) and in the solutions with low initial pH 
(table 1 1); the relatively high values of p/P(OCP) 
in the longest equilibrations indicate that the rate of 
hydrolysis approximated that of dissolution. The 
p/P(OCP) values obtained from the solutions 
which had not reached saturation in the early 
stages of the equilibration and from the solutions in 
which hydrolyzed OCP had increased the value of 
p//'(OCP) were not included in the results shown 
in table 1. The solutions were not considered to be 
in equilibrium with OCP when the p/P(OCP) val- 
ues of the initial solutions or the final solutions 
were 0.2 units higher than the average p^fiTspCOCP) 
value. 



Table 6. Solubility of OCP at 4 °C 









Composition 


of solution 








Time of 




pH 


[Ca] 




[P] 




Ca 


P^,p 


equilib. 






mmol/L 




mmol/L 


P-Po 




(days) 






























0.236 








1 




7.11 


0.459 




0.613 




1.22 


48.37 


4 




7.11 


0.478 




0.618 




1.25 


48.46 


5 




7.10 


0.471 




0.64 




1.17 


48.48 














0.543 








5 




6.77 


0.781 




1.22 




1.12 


48.18 














0.825 








11 




6.63 


0,980 




1.657 




1.18 


48.03 














0.894 








1 




6.52 


1.05 




1.401 




2.07 


48.57 


2 




6.54 


1.05 




1.667 




1.36 


48.29 


4 




6.55 


0.989 




1.664 




1.28 

av 


48,34 
48.26±0,19 


Table 7, 


Solubility of OCP at 4.8 


°C 














Initial 








Composition 


of solution 








conditions 




Time of 
equilib. 


pH 


[Ca] 


[P] 


Ca — Cao 




pH 


HjPO, 


Ca(OH)2 


■?K^ 




mmol/L 


mmol/L 


(days) 




mmol/L 


mmoI/L 


P-Po 




- 


0.248 





11 


6.93 


0.577 


0.791 


1.06 


48.47 




0.543 





1 


6.77 


0.765 


1.23 


1.11 


48.16 




0.849 





17 


6.38 


1.17 


2.05 


0.97 


48.52 




i.7 





17 


6.14 


1.75 


3.21 


1.16 


48.46 


6.08 


1.48 


0.79 


12 


6.44 


1.06 


2.50 


0.26 


48.20 


5.62 


3.0 


1.61 


12 


6.15 


1.95 


3.31 


1.10 
av 


48.21 
48.34+0.16 
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Table 8. Solubility of OCP at 6.0 °C 



Composition of solution 



Time of 


pH 


[Ca] 
mmol/L 


[P] 

mmo!/L 


VK>^ 




equilib. 


Ext. DH 


Davies 





5.92 


4 


3 






5 min 


5.95 










2h 


5.95 


4.19 


3.15 


48.15 


48.19 


Id 


5.97 


4.01 


3.14 


48.13 




4d 


5.95 


3.93 


3.11 


48.26 




7d 


5.97 


3.836 


3.05 


48.30 
av 48.21 ±0.08 


48.34 



Table 9. Solubility of OCP at 18 °C 



Composition of solution 



Time of 
equilib. 

(days) 



pH 



[Ca] 
mmol/L 



[P] 

mmol/L 



P'^sp 



9 

22 



6.76 
6.47 



0.586 
0.840 



0.897 
1.486 



48.24 
48.27 



Table 10. Solubility of OCP at 23.5 °C 





Initial 








Composition 


of solution 








conditions 




Time of 
equilib. 


pH 


[Ca] 


[P] 


Ca-Cao 




pH 


H3PO4 


Ca(OH)2 


pK^ 




mmol/L 


mmol/L 


(days) 




mmol/L 


mmol/L 


P-Po 






0.849 




17 


6.09 


1.62 


2.04 


1.36 


48,31 




1.7 




17 


5.94 


1.74 


3.20 


1.16 


48,37 




3.39 




12 


5.57 


2.94 


5.87 


1,19 


48,67 


6.08 


1.4S 


0.79 


12 


6.26 


0.99 


1.92 


0.45 


48.40 


5.62 


3.0 


1,16 


12 


5.92 


1.69 


3.60 


0.13 

av 


48.40 
48.43+0.14 



Our values of p^sp(OCP) (second column of 
table 1) increased only slightly in the temperature 
range 4 to 37 °C. Those reported by Shyu et al. [5] 
and Heughebaert and Nancollas [6] also varied 
only slightly in the range from 25 to 45 °C, Consid- 
ering the experimental errors and the uncertainty 
of equilibrium constants used, these two studies in- 
dicate that the solubility product of OCP is not 
affected significantly by temperature. The two sets 
of data in the temperature range where they over- 
lap are similar after correction for differences in 
ionic models and equilibrium constants, although 
our values for pATspCOCP) are smaller than theirs. 
An attempt is made in the discussion to compare 
the effect of OCP hydrolysis and the different ionic 



models and equilibrium constants on the 
p/i:,p(OCP) values. 



Discussion 

The OCP is not the most stable phase under the 
conditions studied here or reported in the litera- 
ture. Thermodynamically it will hydrolyze to other 
calcium phosphates; the kinetic rates depend on 
temperature, pH, solid-to-liquid ratio, and calcium- 
to-phosphate ratio in solution. Similar problems of 
instability have also been encountered in the solu- 
bihty studies of carbonate- and fluoride-containing 
apatites [17] and tetracalcium phosphate. The 
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Table 11. 


pjf.p of OCP under Ni 


atmosphere at 37 °C 












Composition of solution 






Time of 


pH 


[Ca] 


[P] 


Ca 


P^.P 


equilib. 




mmo!/L 


mmo!/L 


P-Po 












0.242 






2h 


6.73 


0.331 


0.524 


1.174 


48.93 


2d 


6.45 


0.549 


0.941 


0.785 


48.66 


6d 


6.10 


0.876 


1.64 


0,626 


48.82 
av 48.80 










0.463 






2h 


6.33 


0.633 


1.07 


1.048 


48.80 


6h 


6.28 


0.644 


1.15 


0.967 


48.86 


11 h 


6.14 


0.792 


1.41 


0.840 


48.97 


1 d 


6.03 


0.989 


1.77 


0.758 


48.86 


2d 


5.64 


1.93 


3.71 


0.594 


48.80 
av 48.86 





2.96 





0.935 






2h 


5.81 


1.52 


2.16 


1.241 


48.99 


Id 


5.81 


1.46 


2.16 


1.192 


49.03 


2d 


5.84 


1.49 


2.22 


1,160 


48.83 


5d 


5.82 


1.64 


2.36 


1.151 


48,72 


7d 


5.85 


1.59 


2.32 


1,148 


48.63 
av 48.84 





6.08 


0.79 


1.48 






12 d 


6.26 


0.99 


1.920 


0.45 


48,45 





5.62 


1.61 


3.00 






12 d 


5.70 


1.95 


4.28 


0.27 


48.44 










3.39 






12 d 


5.52 


2.76 


6.02 


1.05 


48.40 










3.50 






2h 


5.37 


2.86 


5,67 


1.32 


49,02' 


Id 


5.42 


2.95 


5,72 


1.33 


48,72 


2d 


5.33 


3.33 


6,47 


1.12 


48,83 


3d 


5,31 


3.41 


6.68 


1.07 


48,86 


7d 


5.24 


4.32 


8.47 


0.87 


48.59 
av 48.75 



' Was not included in solubility calculation. 

longitudinal data given in tables 2-5 and 1 1 clearly 
reveal that hydrolysis was taking place in these 
studies: (1) The pH values of the solutions de- 
creased monotonically with time. This is the result 
of OCP hydrolyzing into a more basic salt. In sev- 
eral instances, the pH changed more than a unit in 
its value without a large change in the p/i'(OCP). 
The lack of change in the value of pArsp(OCP) is 
apparently related to the hydrolysis process being 
slow compared to the rate of dissolution of OCP 
when relatively high solid-to-solution ratios are 
used. (2) The calcium and phosphate concentra- 
tions increased with time after equilibrium had 



been reached. This is in accord with the phase dia- 
gram for OCP. Despite the relatively large changes 
in pH and calcium and phosphate concentrations, 
the p/P{OCP) remained quite constant in most in- 
stances. This is in accord with a process in which 
the composition shifts along the OCP isotherm to- 
ward higher concentrations and lower pH values 
[18] so that the solution remains essentially satu- 
rated. One reason for this behavior is our use of a 
relatively high soUd-to-solution ratio in the equili- 
brations; this favors a high dissolution rate com- 
pared to the hydrolysis rate. 



620 



Volume 93, Number 5, September-October 1988 

Journal of Research of the National Bureau of Standards 



Heughebaert and NancoUas [6] state: "The mean 
molar calcium/phosphate ratios of the solid phases 
after equilibrium in these experiments were 
1.30+0.05 at 25 °C and 1.33±0.02 at 45 "C con- 
firming the good stabihty of OCP in the aqueous 
media under these conditions of temperature and 
pH. In particular, hydrolysis of OCP into an HAp- 
like phase at 45 °C in the pH range 5.5—6.8 was 
never detected under the experimental conditions 
used." The changes in the solution composition at 
25 °C [6], which is much more sensitive to the hy- 
drolysis as compared to the changes in Ca/P ratios 
of the solid, indicate some hydrolysis of OCP. The 
ACa/AP values calculated from the Heughebaert 
and NancoUas 25 °C data (table II, ref. [6]) reveal 
that the ACa/AP ratios were larger than 1.33 (table 
12) when equiUbrium was approached from super- 
saturation (av ACa/AP =1.53 when the unreUable 
value from expt. 52G is omitted) and all the ACa/ 
AP ratios were smaller than 1.33 when equilibrium 
was approached from undersaturation. There is no 
overlapping in the two sets of data. The above re- 
sults are in accord with hydrolysis processes in 
which a product more basic than OCP is being 
formed. For example, consider a process in which 
precipitation of OCP alone is occurring. The cal- 
cium lost from solution divided by phosphate that 
is lost in the same time interval, ACa/AP, would be 
1.33. Now if, say, OHAp precipitated along with 
the OCP (or was formed by hydrolysis of OCP in 
situ), the quantity ACa/AP would be larger than 
1.33. This is the situation with most of the 25 °C 
data from supersaturation (table 12). The reverse 
would be the case when approaching equilibrium 
from undersaturation. If only OCP were dissolv- 

Table 12. Calculated ACa, AP, and ACa/AP values from the 
25 °C data of Heughebaert and NancoUas [6]. The equilibrations 
44G through 52G are equilibrations approached from supersatu- 
ration; those labeled 44D to 52D are from undersaturation 



# 


ACa 


AP 


ACa/AP 


52G 


0.014 


0.005 


2.8 


50G 


0.20 


0.016 


1.25 


48G 


0.37 


0.22 


1.68 


47G 


0.26 


0.20 


1.30 


46G 


0.27 


0.15 


1.80 


44G 


0.18 


0.11 


1.64 


52D 


0.092 


0.096 


0.96 


SOD 


0.157 


0.185 


0.85 


48D 


0.151 


0.187 


0.81 


47D 


0.197 


0.206 


0.96 


46D 


0.215 


0.227 


0.95 


44D 


0.266 


0.236 


1.13 



ing, ACa/AP would be 1.33. But if some OHAp 
were forming, Ca and P would be removed in a 
ratio greater than 1.33, leaving a net dissolved 
ACa/AP less than 1.33. Thus, the ACa/AP data at 
25 °C (table 12) indicate that hydrolysis of OCP 
had occurred. The ACa/AP ratios calculated from 
the data at 37 °C by Shyu et al. [5] also indicate the 
occurrence of OCP hydrolysis. On the other hand, 
the solution compositions of the 45 °C equilibra- 
tions [6] do not indicate that hydrolysis of OCP 
had occurred, and the TpKsp(OCP), calculated from 
these data using our ionic model, is 48.6±0.4, the 
same as our value at 37 °C. 

The pJsrsp(OCP) values for four groups of investi- 
gators using different ionic models and equilibrium 
constants are compared in table 1. It is apparent 
that the p^5p(OCP) values reported by Heughe- 
baert and NancoUas [6] and by Shyu et al. [5] at 25 
and 37 °C are about an order of magnitude higher 
than those reported here. Our 37 °C value is higher 
than the values attributed to Madsen [7]. Madsen's 
two values for p^sp(OCP) are recalculations from 
his calculated solubility data that were based on an 
equilibrium model which did not mclude the ion 
pairs CaH2P04+, CaHPO*'', and CaOH". The 
Moreno et al. [8] value at 25 °C is based on a single 
composition (the singular point for brushite and 
OCP). Thus, even though it has been recalculated 
to take into account calcium and phosphate ion 
pairs, relatively little reliance can be placed on this 
value [8]. We discuss the causes for the discrepan- 
cies between our p^sp(OCP) values and those of 
Shyu et al. [5] and of Heughebaert and NancoUas 
[6] in the following sections. 

Thermodynamic Systems and Ionic Model 

Heughebaert and NancoUas [6] and Shyu et al. 
[5] used a five component system, Ca(OH)2 
-H3PO4-KNO3-H2O-KOH. The component 
KNO3 was included in the Heughebaert and Nan- 
coUas [6] and Shyu et al. [5] systems so that extra- 
polations could be made to zero ionic strength in 
evaluating pK^piOCP). The component KOH was 
used to adjust the pH. We used mostly the ternary 
system, Ca(OH)2-H3P04-H20, except for the 
studies with initially supersaturated solutions. Our 
use of a ternary system was motivated by the desire 
to keep the system as simple as possible so as to 
avoid unanticipated ionic interactions. In the 
ternary system, the ionic strengths were generally 
so low that the ion activities could be calculated 
directly using the extended Debye-Hiickel 110] or 
Davies equation [20]. However, both methods 
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should yield satisfactory results provided they are 
based on a valid ionic model of the system. It is 
interesting that the p^fiTspCOCP) values are not signif- 
icantly different when either the extended Debye- 
Huckel equation or Davies equation is used as 
shown in tables 1 and 8. 

Equilibration Conditions 

Heughebaert and Nancollas [6] equilibrated their 
samples for 25 days at 25 °C and 12 days at 45 °C, 
apparently without sampling at shorter time peri- 
ods; Shyu et al. [5j sampled at termination of each 
experiment after periods of 4, 10, 30 or 40 days. As 
indicated in the Methods section, we used higher 
solid-to-solution ratios (10 mg/mL) as compared to 
Shyu et al. [5] and Heughebaert and Nancollas 
(3 mg/mL), and sampled frequently at much 
shorter time periods, In doing so, a constant value 
of pArsp(OCP) was obtained in a relatively short 
time. In the experiments with prolonged equilibra- 
tion times, a subsequent increase in the value of the 
apparent pArsp(OCP) was observed which we at- 
tribute to hydrolysis of the OCP to a more basic 
and less soluble form of calcium phosphate. Table 2 
clearly shows increases in the apparent p^sp(OCP) 
after the longer equilibration times. 

Heughebaert and Nancollas [6] and Shyu et al. 
[5] approached equilibrium from both undersatura- 
tion and supersaturation. Most of our equilibrations 
were from undersaturation, but good agreement 
exists between our data from supersaturation and 
undersaturation. 

The value of a pK^(OCF) calculated from a set 
of experimental data is particularly sensitive to 
small errors in the measurement of pH. Glass elec- 
trodes were used by both groups of investigators. 
Heughebaert and Nancollas [6] and Shyu et al. [5] 
used a silver/silver chloride reference electrode 
which incorporated an intermediate liquid junction 
containing potassium chloride solution at the same 
ionic strength as the solution being studied and. 



therfore, avoiding the diffusion of the potassium 
chloride into solution. We used a standard calomel 
reference electrode and minimized the diffusion of 
the potassium chloride by avoiding prolonged con- 
tact between the electrode and the solution. 

Equilibrium Constants 

The p^sp(OCP) is calculated from the pH, {Ca], 
and [PO4] data through the use of an ionic equi- 
librium model [15]. This model uses dissociation 
constants of H3PO4 and H2O and the association 
constants of the ion pairs CaHP04, CaHiPO^, and 
Ca(OH)+ and, in this study, CaCOj and CaHCOj. 
A comparison of the constants used by Heughe- 
baert and Nancollas [6] with those used by us 
(shown in table 13) reveals substantial differences 
in the formation constants for the ion pairs at some 
of the temperatures. In table 1 are given the 25, 37 
and 45 °C values of p^sp(OCP) reported by Shyu et 
al. and Heughebaert and Nancollas [5,6]. Below 
each of these values is given in parentheses the 
value we calculated using their data for pH, [Ca], 
[PO4], and neutral ions with our model and equi- 
librium constants. The pA'sp(OCP) values calcu- 
lated by Heughebaert and Nancollas [6] are 
significantly higher than our values calculated 
from the same experimental data; these differences 
are probably due to the different equilibrium con- 
stants used, since both the extended Debye-Hiickel 
equation and the Davies equation gave the same 
results. Clearly, there is a need to reassess the val- 
ues of the ion pair formation constants since they 
contribute significantly to the differences in the 
pJirsp(OCP) values as shown in table 1. 

The pArsp(OCP) values in our study are similar to 
those calculated from their data [5,6] with our 
model and equilibrium constants (table 1); our val- 
ues are slightly lower than theirs at 23.5 to 37 °C 
probably due to the hydrolysis of OCP in their 
study. 



Table 13. Equilibrium constants [5,6,15] 



25 'C 



37 "C 



45 °C 



This study Heughebaert et al. [6] This study Shyu et al. [i] 



HjPO4:^H++H2P0r 
HiP04:;i;H++HP042- 
HP04^-:?±H++P04^- 
Ca^+ + HiP04-;?=CaHjP04+ 
Ca^+ +HP04^-=:CaHP04 
Ca^+ + P04^"^=iCaP04- 

Cu^'-.|.oifcCuori+ 

H++OH-=iHjO 



This study 


Heughebaert et al, [6^ 


5.63X10-' 


5.63x10-5 


6.60x10-* 


6.61X15-= 


9.02x10-'^ 


6,61X10-" 


5.57 


36.5 


503 


787 


2.9X10' 


3.86X10' 


20,0 


26.4 


4.04XI0-" 


4.09X10-" 



7.11x10-5 

6.31xlO-» 

4.52X10-" 

S.4S 

264 

2.9x 10* 

20.0 

1.013X10- It 



7.11X10-' 

6.30X10-' 

4.73 X 10-" 

25.6 

548 

2.9X 10' 

13,8 

1.004X10-" 



6.22X10-' 6.22x10-5 



6.58x!0-» 

6.84X10- '5 

7.01 

355 

2.9X10' 

20,0 



6.58X10-' 
6.61x10-" 

31.9 
631 

3.46x10' 
21.3 



2.41x10-'* 2.42X10-'* 
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Appendix 

The ionic activity coefficients were calculated 
from the extension of the Debye-Hiickel limiting 
law [19] or Davies equation [20]. The expanded 
Debye-Hiickel equation is 

log^= -^z V//(l +£fly/) 

and the Davies equation is 

logy?= -^z=[V//(i + V/)-o.3/], 

where f\ and / are the activity coefficients of spe- 
cies i and the ionic strength, respectively. A and B 
are functions of the temperature and the dielectric 
constant of the solvent (Appendix table 4, ref. [19]). 
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For the parameters Cj we used for Ca^"*": 6x10 ^ 
cm; for H2PO4-: 4.5x10-* cm; for HPOr and 
PO4-: 4X10-' cm; for H+: 9x10-' cm; and for 
OH-; 3.5X10-^ cm [21], The use of different val- 
ues of Cj thermodynamically violates the Gibbs- 
Duhem equation [22]. The equilibrium constants 
for ion pairs, phosphoric acid and water are shown 
in table Al [15,16,22]: logiJ<:= A^/T^+Az/ 
T+A^+A^T+AslogioT, where r=kelvins. The 
temperature dependence of the k^ value 
for H3PO4 is given by: — log ioA = 12.45 
-0.015(°C-18.0), where °C is the temperature. 

Table Al, Parameters for equilibrium constants 



Ion pair 


^1 


^1 


^3 


A, 


As 


H3PO4 kt 


0. 


-799.31 


4.5535 


-0.01349 


0. 


H3PO, k. 


0. 


- 1989.32 


5.4374 


-0.02001 


0. 


K 


0. 


-4470,99 


6.0875 


-0.01706 


0. 


H2CO3 k, 


-1684915. 


21834.37 


-356.3094 


-0.06092 


126.8339 


HjCOj ki 


-563714. 


5151.79 


-107.8871 


-0.03253 


38.9256 


H2CO3 k„ 


669365. 


-6919.53 


108.3865 


0,01985 


-40.4515 


(gas partition) 












CaH2p04 


0. 


-8232.51 


57.1388 


-0.09592 


0. 


CaHPO, 


0. 


12213.79 


-81.1027 


0.14275 


0. 


CaPOi 


0. 


0. 


6.4624 


0. 


0. 


CaOH 


0. 


0. 


1.3010 


0. 


0. 


CaHCOj 


0. 


-34765.05 


1209.1200 


0.31294 


478.782 


CaCOj 


0. 


35512.75 


-1228.7320 


-0.29944 


-485,818 
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